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The Gluon Propagator in Momentum Space 
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We consider quenched QCD on a 16 3 x 40 lattice at j3 = 6.0 . We give preliminary numerical results for the 
lattice gluon propagator evaluated both in coordinate and momentum space. Our findings are compared with 
earlier results in the literature at zero momentum. In addition, by considering nonzero momenta we attempt to 
extract the form of the propagator and compare it to continuum predictions formulated by Gribov and others. 



1. INTRODUCTION 

The possibility of studying nonperturbatively 
on the lattice gauge-dependent quantities pro- 
vides in principle a unique tool to test QCD at 
the level of the basic fields entering the contin- 
uum Lagrangian. From this point of view, the 
gluon propagator in the quenched approximation 
is perhaps the simplest quantity. From its study 
one expects to obtain among other things a bet- 
ter understanding of the infrared behaviour of the 
theory and of the mechanism of gluon confine- 
ment. 

The nonperturbative behaviour of the Eu- 
clidean gluon propagator has been investigated 
in the continuum by many authors with different 
methods and in different gauges |l|-|^|. In partic- 
ular, a very peculiar momentum dependence has 
been predicted as a consequence of a modifica- 
tion of the standard path integral Faddeev-Popov 
formula in the Landau gauge by the introduction 
of a nonperturbatively correct gauge-fixing proce- 
dure |l],Q. Such improved implementation of the 
Landau gauge is expressed by the equations 



where FPh4] is the Faddeev-Popov operator in 
the Landau gauge, which in general is not positive 
definite. The positivity requirement in ([I]) can be 
seen as a recipe to get rid (although not com- 
pletely, see for example of Gribov copies 0. 
In the gauge (Q) , the (transverse) gluon propaga- 
tor in momentum space has been argued to be of 
the formEir 



(2) 



^Gribov (k) 



k 4 + b 4 



where & is a dynamically generated mass scale. 
The form (Q) for the propagator in momentum 
space implies that in the continuum 



GGribov(k — 0, t) 
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FP[A] > (1) 



(3) 
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Remarkably, the same predictions were also ob- 
tained in the study of Schwinger-Dyson equa- 
tions (|]. 

2. THE LATTICE PROPAGATOR 

The lattice gluon field can be defined asfjj: 

U,(n)~Ul(n) I ( U,(n)-Ul(n) \ 
2fa 3* 2fa (4) 
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where a is the lattice spacing. Thus the lattice 
gluon propagator in space is the expectation 
value of: 



G^(x,y) = Tr(A^(x) A v {y)) 



(5) 



An important point is that on the lattice one 
can define and implement the analogue of the 
gauge condition (|l|) and reobtain from analyti- 
cal arguments the predictions for the propagator 
mentioned in the above section[|). 

In fact, given any link configuration {[/}, one 
can define a function of the gauge transformations 
g on {U} 

F u{9} = ~J2 Re Tr (W + U ^ n ~ A)) . ( 6 ) 



where V is the lattice volume and U 9 in- 
dicates the gauge-transformed link U 9 (n) = 
g(n)U IM (n)g^ (n + p.). An iterative minimization 
of Fjj[g] obtained by performing suitable gauge 
transformations generates a configuration {U 9 } 
that satisfies the lattice version of (|), defined m 
terms of a lattice Faddeev-Popov operator. Then 
it is natural to try and test numerically predic- 
tions like (||) and (||). 

Numerical studies have been performed in 
the past years for the zero spatial momentum 
Fourier transform of (ff), namely G(k = 0,t) = 
Si=i Gu (k = 0, t) These studies reported 
some evidence of an effective gluon mass that in- 
creases with the time separation. This feature, 
which would be unacceptable for the propaga- 
tor of a real physical particle since it violates the 
Kallen-Lehmann representation, is in qualitative 
agreement with the continuum prediction ^) and 
may be in principle acceptable for a confined par- 
ticle§,§. 

Our work aims to test at a more quantitative 
level continuum predictions and to extend the 
above results through the study of the gluon prop- 
agator at nonzero momenta. By requiring consis- 
tency between zero and nonzero momentum re- 
sults, one has a better chance to determine the 
propagator's analytical form. 

3. NUMERICAL RESULTS 



3.1. Technical Remarks 

It is perhaps worth remarking that, unlike sim- 
ulations involving quenched quark propagators, 
evaluations of purely gluonic correlation functions 
can take full advantage of the translational sym- 
metry of the theory in order to improve statistics. 
On the other hand, such quantities turn out to be 
very sensitive to the numerical accuracy of gauge 
fixing. Empirically, it turns out that when the 
minimization of -FVlff] has reached an accuracy 
w .05% the signal for the propagators is suffi- 
ciently stable against additional gauge fixing. 

3.2. Results 

We report results for a set of 25 configurations 
on a 16 3 x 40 lattice at (3 = 6.0. As a first step 
we have evaluated G(k = 0,t); our results con- 
firm that the propagator exhibits a massive de- 
cay in time with an effective mass a * m(t) = 
ln( G i fc ~°' t ) ) that increases with t. In Fig.Rl) we 

plot a * m(t) versus t in lattice units with jack- 
knife errors. Assuming the value of the inverse 



Figure 1. Effective gluon mass in lattice units 



lattice spacing a -1 w 2.0GeV at f3 = 6.0, the 
effective gluon mass that we measure ranges ap- 
proximately between 200 and 800 MeV. 

We first attempt a 2-parameter least-squares 
fit of our data to the continuum form (|^) without 
taking into account the correlations in the data. 
The parameters are an overall normalization fac- 
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tor and the mass scale b. The results are given in 
Fig. H and show a very good agreement between 
the data and the fitting points. For a -1 = 2 
GeV one obtains b — 225 ± 5 MeV, where the 
quoted error comes from the covariance matrix of 
the fitting parameters and does not include the 
systematic error on a -1 . 

Good agreement is also obtained by using 
the form commonly referred to as particle + 
ghost, that is G(k = 0,t) » C x exp{-Mit) + 
C 2 exp(— M 2 t), where C 2 is constrained to be neg- 
ative. 

On the other hand, one cannot get good agree- 
ment if one uses a conventional 4-parameter dou- 
ble exponential form, that is if one constrains C 2 
in the above formula to take positive values. In- 
deed, in this case the effective mass would always 
decrease with t, in contrast to what is observed. 



Figure 2. G(k = 0,t) (error bars shown) and fit 
to the form Gcribov (diamond points) 



It is well known that the data points obtained 
from a Monte Carlo simulation are in general 
statistically correlated; in the present case, the 
correlated data are the values of the propagator 
G(k = 0, t) at different timeslices. 

The standard way to take into account this ef- 
fect when performing \ 2 fits is to use the defini- 
tion of x 2 that involves the full covariance ma- 
trixjn]]. Such a definition reduces to the stan- 
dard one when the covariance matrix is diagonal, 



which would happen if the data points were un- 
correlated. 

By inspection of the covariance matrix for 
G(k = 0, t) it turns out that the off-diagonal ma- 
trix elements are typically of the same size as the 
diagonal ones, i.e. our data points are highly cor- 
related in t. Consequently, when we perform \ 2 
fits taking into account the full covariance ma- 
trix, the fits are not well controlled because the 
correlation matrix is nearly singular. However, 
we still find that Gcribovik = 0,t) fits the data 
better than other forms. There is also qualitative 
agreement between our results for G{k = 0, t) and 
previous ones from other groups 

In spite of the difficulties in the statistical anal- 
ysis, our interpretation of the results for G(k = 
0, i) receives a strong support from the analy- 
sis of the momentum space propagator G(k) = 
X^=i G^(k). It turns out that such a quantity 
is very well determined in a significant interval 
of physical momenta, ranging from the lattice in- 
frared cutoff k a — to k 3k Q (see Fig. ||). 
In this range we fit the data to the continuum 
formula @ and, for a comparison, to a standard 
massive propagator G mass {k) ~ k z+ m s ■ 

An interesting point is that the covariance ma- 
trix associated with G(k) turns out to be much 
more "diagonal" than the one for G(k = 0,t); in 
other words, the data points are much less corre- 
lated in momentum space than they are in t. As 
a consequence we have been able to obtain good 
fits for G(k) with or without taking into account 
correlations. We show in Fig. | a fit of G(k) 
to the form Gcribovik). With the full covariance 
matrix, we get for the fit in Fig. Hxi ff r = 1-5 aim 
b = 322 ± 8 MeV, assuming again a" 1 = 2 GeV 
and neglecting systematic errors. 

We compare this result to the best fit that one 
can obtain from the standard massive propaga- 
tor, for which we obtain Xd gr — 2-9. On the other 
hand, the b values that one obtains from the fits 
in coordinate and momentum space differ signif- 
icantly, since they are respectively b = 225 ± 5 
and b — 322 ±8 MeV. Since an appreciable differ- 
ence also occurs for the mass parameter m when 
we compare momentum and x— space fits to a 
simple massive propagator, we think that such 
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requires the study of the propagator at very low 
momenta, in order to observe the suppression pre- 
dicted by (^|). Of course such a study calls for very 
big lattices. 

The work in progress aims to obtain first a bet- 
ter understanding of systematic and statistical er- 
rors. After that, a study of the scaling properties 
of the mass scale associated to the gluon propa- 
gator is in order and, in a later stage, the issue of 
the gauge dependence of the propagator will be 
addressed. 



Figure 3. momentum space propagator vs. 
GeV (assumes 1/a = 2.0 GeV) 



Figure 4. G(k) (error bars shown) and fit to the 
form Gcribov (diamond points) 



discrepancies may be related to the different role 
that finite-size effects play in the two calculations. 
Further investigation of this issue is in progress. 

4. CONCLUSIONS 

We think that our results provide a significant 
(although preliminary) check of the continuum 
predictions (||) and (||). In particular, the study 
of the propagator in momentum space appears 
very promising since the data for such quantity 
turn out to be statistically rather clean. 

Recalling M) it is clear that a conclusive test 
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